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I 

, The triaxial deformation in ^''Be is investigated using a microscopic 

\ a+a+n+n model. The states of two valence neutrons are classified based 

on the molecular-orbit (MO) model, and the vr-orbit is introduced about 
the axis connecting the two a-clusters for the description of the rotational 
. bands. There appear two rotational bands comprised mainly of K"^ = 0"*" 

ff^ , and K'^ = 2+, respectively, at low excitation energy, where the two valence 

\ neutrons occupy K'^ = 3/2" or K'^ = 1/2" orbits. The triaxiality and the K- 

' mixing are discussed in connection to the molecular structure, particularly, to 

the spin-orbit splitting. The extent of the triaxial deformation is evaluated in 
terms of the electro-magnetic transition matrix elements (Davydov-Filippov 
model, (^'invariant model), and density distribution in the intrinsic frame. 
The obtained values turned out to be 7 = 15° ~ 20°. 
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I. INTRODUCTION 



Recently, exotic structures of the Be isotopes have been theoretically and experimentally 
studied, and many new phenomena have been discussed One of them is the appearance 

of the cluster rotational band structure in the excited states of ^°Be (a+^He) 0,0] and ^^Be 
(a-l-^He, '^He+'^He) M. We have suggested [§-0], based on the molecular-orbit (MO) 



model, that the development of a-a cluster structure depends on the neutron orbits located 
around the core comprised of a-clusters. In the second 0"^ state of ^°Be, an anomalously 
prolonged a-a clustering structure emerges due to the valence neutrons located along the 
a-a axis 

The nucleus ^°Be has been known to have a strong (3 deformation due to the a-a core. 
In this paper, we discuss another aspect, a triaxial deformation. The triaxial deformation 
is possible as a result of the dynamics of the two valence neutrons. The triaxiality of ^°Be 
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has been theoretically discussed based on the deformed oscillator model [|TT|, in which a 7 



distortion of 34.8° has been predicted. According to Davydov-Filippov model [|12|, this 7 
value suggests that the excitation energy of the second 2"*" state is 2.3 times higher than that 
of the first 2"*" state. In ^°Be, the first 2"*" state is observed at 3.358 MeV, and several 2"*" 
states have been observed around 8 MeV region. Because of the presence of the 2^ states in 
this energy region, recently, an experimental study has been performed in order to identify 
a triaxial structure of ^"Be [P^. 

In this paper, we discuss the triaxiality of ^°Be as a composite system of the a-a core 
and two valence neutrons, comparing with other theoretical models. In our previous study 
on ^°Be 0, all of the observed low-lying positive- and negative-parity states are explained 
as combinations of three basic orbits [K'^ = 3/2~, 1/2+ , and 1/2") of two valence neutrons 
around the two a-clusters. Here, the z-axis is taken to be the axis connecting two a-clusters. 
If we adopt K'" = 3/2~ or l/2~ orbits for the two valence neutrons, there appear two 
rotational bands in low energy; one is dominated by the K = intrinsic structure and the 
other hj K = 2. The calculated two 2+ states of these bands can be related to the observed 
first 2"*" state at 3.358 MeV and the second 2+ state at 5.958 MeV. It is therefore important 
to show how the triaxial intrinsic configuration emerges in these states, and how the orbits 
of the valence neutrons deviate from the axial symmetry. Here, we calculate the electro- 
magnetic transition between these K = and K = 2 bands (B(E2: K = 2 ^ K = 0)) 
as a signal of a triaxial deformation. This transition is suppressed when the orbitals are of 
pure axial symmetry. However, the orbitals may deviate from the axial symmetry, when the 
valence neutrons are mutually more correlated and form a localized di-neutron pair due to 
the neutron-neutron interaction. The recoil effect of the valence neutrons with respect to 
the a-a core then plays a role to break the axial symmetry of the charge distribution. 

This paper is organized as follows. In section II, we give a description of the single- 
particle orbits around the two a-clusters based on MO. In section HI, we show our results 
on ^°Be, and a conclusion is given in section IV. 



II. EXTENDED MOLECULAR-ORBIT MODEL 
A. Wave function of ^'^Be 
We introduce a microscopic a+a+2n model for ^°Be. The neutron configurations are 



introduced based on the molecular orbit (MO) picture |14-16|. The total wave function is 
fully antisymmetrized and expressed as a superposition of Slater determinants with various 
configurations of the valence neutrons. The Slater determinants are also superposed with 
respect to different relative distances between the two a clusters {d). The projection to the 
eigen states of angular momentum J is numerically performed. All nucleons are described 
by Gaussians with a common oscillator parameter s = set equal to 1.46 fm. The a 
cluster located at Ra on the 2;-axis contains four nucleons. 

Here, G represents a Gaussian: 

3 

Gr. = (^yexp[-u{r-R^)% u=l/2s', (2) 
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where Ra = Ra^z with being the unit vector pointing the z-direction. The actual value 
of Ra is d/2 or —d/2 in this work. The wave function of the i—th valence neutron {(piXi) is 
expressed by a linear combination of local Gaussians: 

(t)iXi = ^9jGRjXi-, (3) 



2z/\3 



Gr, = [-) exp[-Krl-i?,n, (4) 

where Rj is a parameter corresponding to the Gaussian center. 

The two valence neutrons are classified according to the MO picture. In the MO model, 
the states of the valence neutrons are expressed by a linear combination of orbits around 
two a clusters. Here, we notice that the valence neutrons and the neutrons in the a clusters 
are identical particles, and the valence neutrons must be orthogonal to the forbidden space 
due to the a-a core. The lowest orbit then has one node and negative parity, that is the 
p-orbit. We take the harmonic-oscillator- type wave function for the p-orbit. Such p-orbits 
in the x-, y-, and z-direction are denoted as ipx, ipy, and ipz, respectively. We now classify 
single particle orbits for the valence neutrons in terms of K quantum numbers: 

ipx + ii^y = {x + iy) exp[— z/r^] oc rYu exp[— z/r^], (5) 
ipz = z exp[—i'r'^] Qc rYio exp[— i^r^], (6) 

i'x — i^y = {x — iy) exp[— z/r^] oc rFi_i exp[— z/r^]. (7) 

These p-orbits have K = 1, 0, and —1, respectively. Here, x, y, and z are relative to Ra- 

Now we construct MO from these p-orbits. Since each valence neutron can move around 
one of the two a-clusters, there should be two sets of the p-orbits defined with respect to 
those two a-clusters. Thus, the {ipx ± ii^y) orbit whose center is shifted at +a (—a) on the 
z-axis is denoted as {ipx ± ii^y)+a {{'ipx ± i4'y)-a)- As linear combinations of these orbits, the 
lowest MOs are expressed as: 

^ = i^x + i^Py)+a + i^x + i^y)-a, (8) 
= {^^ - ^A)+a + (V^x - ^^y)-a. (9) 

These orbits clearly have K = 1 and K = —1, respectively. In these cases, the classical 
picture of p-orbit is a circular motion about the a-a {z) axis. This is the so-called vr-orbit. 

If we take a linear combination of {%l)z)+a and {il)z)-ai the orbit becomes the so-called a- 
orbit, just along the a-a [z) axis. This orbit is a higher-nordal orbit and only relevant to the 
second O"*" state, discussed also in antisymmetrized molecular dynamics (AMD) calculations 
r?! and in stochastic variational method (SVM) calculations [jl8[, and now we disregard 



this orbit. 
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These {iljx)+a, {4'x)-a, i4'y)+a, and {4'y)-a orbits can be approximated by a combination 
of two local Gaussians, whose centers are shifted by a variational parameter b perpendicular 
to the z-axis. We thus use the following wave functions, 

{4'x)+a OC Gaez+bSx ~ Gaez-be^-, O^y) +a OC Gae^+bSy ~ Gaez-bSyi ■ ■ ■ (10) 

The values of these parameters a and h are variationally determined by using the cooling 
method in AMD |ll9|-p^ independently for each a-a distance. Since the rotational symmetry 



about the 2;-axis is broken with Eq. (10), the wave functions do not have exactly conserved 
quantum number K. On the other hand, the parameter h is small enough usually, and the 
K"^ number is preserved to a good extent, and is expressed hereafter as . 

When these orbital part of the wave functions are coupled with the spin part, the orbits 
with K'^ = 3/2, 1/2, —1/2, and —3/2 are introduced, as follows: 

I3/2-) = T), I1/2-) = i), (11) 

and their time reversal, 

\-3/2-)=^'_{0\nl), I - 1/2-) = T). (12) 



B. Hamiltonian 

The Hamiltonian consists of the kinetic energy term, the central two-body interaction 
term, the two-body spin-orbit interaction term, and the Coulomb interaction term: 

^ = E^^-r- + E^^. + E^.? + E^^(i-^:)(i-^i)- (13) 

i i<j i<j i<j *J 



The effective nucleon-nucleon interactions are Volkov No. 2 |24] for the central part and the 



G3RS spin-orbit term for the spin-orbit part, as follows: 



V,j = {^le-'^^"'^ + V2e-^^''^^}{W - MP^P^ + fiP" - HP^}, (14) 

Vi; = V^'{e~''''l - e-''< }L ■ SP^u (15) 

where P31 is a projection operator onto the triplet odd state, and L and S operators represent 
the relative angular momentum and total spin of the interacting two nucleons, respectively. 
The parameters are Vi = -60.650 MeV, V2 = 61.140 MeV, Oi = 0.309 fm'^ and 02 = 0.980 
fm~^ for the central interaction, and V^'** = 2000 MeV, ai = 5.00 fm~^, and a2 = 2.778 
fm~^ for the spin-orbit interaction. We employ the Majorana exchange parameter M = 0.6 
{W = 0.4), the Bartlett exchange parameter B = 0.125 and the Heisenberg exchange 
parameter H = 0.125 for the Volkov interaction (using B and H, there is no neutron- 
neutron bound state). All these parameters are determined from a + n and a + a scattering 



phase shifts and the binding energy of the deuteron [EQ] 
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C. Configurations of tlie valence neutrons 

For i°Be, we introduce three configurations of $((3/2-)2), $((l/2-)2), and $(3/2" 1/2") 
for the two valence neutrons, and mixing amphtude of these configurations are variationally 
determined after the angular momentum projection. 

It will be shown that the $((3/2~)^) configuration, where valence neutrons occupy orbits 
with K'^ = 3/2^ and K'^ = —3/2^ is the main component of the ground state. The energy 
shift due to the spin-orbit interaction acting effectively on a valence neutron is proportional 
to —{I ■ s). It is, therefore, favored that the orbital part and the spin part of K are coupled 
to be parallel. This means that the spin-up valence neutron {\n t)) has K'^ — 3/2~, while 
the spin-down valence neutron {\n I)) — — 3/2~, giving rise to the wave function of the 
two valence neutrons: 

$((3/2-)^) = ^[#V?^(</'ciXci)(toc2)], (16) 

with 

0ciXci = V'f''|riT), toc2 = V'Tl^i)- (17) 

We also introduce a basis state $((1/2^)^), which the spin-orbit interaction does not 
favor. They are defined by flipping the spin-part of the 3/2^ and —3/2^ wave functions, 
and the valence neutrons occupy orbits of K'^ — 1/2' and — — 1/2~: 

^clXcl^V^f'lni), 0c2Xc2=V'Tl^T)- (18) 

The $(3/2~ 1/2") configuration with K — 2 \s constructed as a combination of the 
valence neutrons in the — ?>/2~ and — 1/2" orbits. 

0ciXci = <''|nT), </>c2Xc2 = <''|ni). (19) 



III. TRIAXIAL STRUCTURE IN ^"BE 

In ^°Be, we adopt three configurations $((3/2^)^), $((1/2^)^), and $(3/2" 1/2"), where 
single particle orbits with ^ = 3/2", 1/2^ have been introduced. The values of the param- 
eters for the valence neutrons are obtained variationally. The parameter a describes the 
positions of the Gaussian centers on the a-a {z) axis, and parameter h corresponds to the 
rotation radius of the 7r-orbit about the a-a axis. The optimized values are listed in Table 
I. 



Table I 



The parameter a is obtained to be approximately the same value as d/2. Therefore, the 
picture "orbit around the a-cluster" works well, while there is a small deviation of the orbit 
along the ^-axis. The parameter h depends on the configurations. For $((3/2~)^), h appears 
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to be '--'I fm. For $((l/2~)^), b becomes longer as being 2 fm, since the spin-orbit interaction 
acts attractively for $((3/2^)^) whereas repulsively for $((1/2^)^). 

Before performing the angular momentum projection, we discuss the triaxial deformation 
of the intrinsic wave function. The intrinsic density of $((3/2~)^) is shown in Fig. 1 ((a): 
x2;-plane, (b): xy- plane), where the a-a distance is chosen to be optimal, as it turns out to 
be 3 fm. 



Fig. 1 (a), (b) 



The intrinsic density is defined as a snap-shot of a rotating system. The plane in which the 
spin-up valence neutrons stays at a given moment is defined as the a;2;-plane, with the 2;-axis 
being the a-a axis. Therefore, the spin-up valence neutron is fixed on the x2;-plane, and the 
spin-down valence neutron occupy the normal K = — 3/2~ orbit. The intrinsic density on 
the xz-pla.iae (Fig.l (a)) shows a large /3 deformation along the z-a.xis due to the presence of 
two a-clusters. As for the xy-plane (Fig. 1 (b)), the density shows deviation from circular 
distribution and mixture of a triaxial component is evident. The contour-line of 0.01 in Fig. 
1 (a) and (b) suggests the deformation of this nucleus is /3 ~ 0.5 and 7 = 11". 

We then perform the angular momentum projection, where the wave function is rotated 
and integrated over the Euler angle, and the rotational symmetry is restored. Here, the 
X-mixing between K — and K — 2 finally determines the degree of the triaxiality of the 
system. 

This intrinsic wave function is numerically projected to the eigen states of angular mo- 
mentum, and the basis states with parameters listed in Table I are superposed by generator 
coordinate method (GCM). Here, the coefficients representing linear combinations of Gaus- 
sians for each single particle orbit are treated as variational parameters to take into account 
deviations from the original orbits. Because of this, not only "single-particle" MO state 
where the valence neutrons independently rotate around the a-a core, but also more com- 
plex states where the valence neutrons are mutually more correlated as a di-neutron cluster 
can be included. 

The energy levels of ^°Be are shown in Fig. 2. Using our framework, the binding energy 
of one a-cluster is calculated to be 27.5 MeV, and the threshold energy of free a-\-a-\-n-\-n 
system is 2 x (—27.5) = —55.0 MeV. Experimentally, the ground state is lower than this 
energy by 8.4 MeV. 



Fig. 2 



The left and middle columns in Fig. 2 show the calculated energy levels with the constraints 
of K = and K = 2, respectively, and the right column shows the result after the fT-mixing. 
In K = column, it can be seen that the ground 0+ state at —60.5 MeV, the 2+ state at 
—56.9 MeV, and the 4+ state at —47.6 MeV form a rotational band structure. They fit 
quite well into the J{J-\-l) rule, and the dominant component is (3/2~)^ for the two valence 
neutrons. For K = 2 column, the 2+ state at —54.9 MeV, the 3"'" state at —51.1 MeV, and 
the 4+ state at —46.1 MeV form a rotational band structure. They also fit quite well into 
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the J(J + 1) rule. The K = 2 band dominantly has a component of (3/2~)(l/2~) for the 
two valence neutrons, where one of them feels the spin-orbit interaction attractively and the 
other feels repulsively. After the mixing of these two bands, the 2"*" state with K = at 
-56.9 MeV becomes -57.3 MeV {E^ = 3.1 MeV, experimentally 3.358 MeV), and 2+ state 
with K = 2is slightly pushed up to -54.9 MeV {E^ = 5.6 MeV, experimentally 5.958 MeV). 
The first 2^ state has the squared overlap with the K = state by 0.96, and the second 2"*" 
state has the squared overlap with the K = 2 state by 0.92. Therefore, the second 2+ state 
has the component of = by 8%, and it can be considered as an indication of a triaxial 
deformation. 

This triaxiality of the 2+ states is reflected in the electro-magnetic transition rate, and 
B(E2) values are summarized in Table II. 



Table II 



The E2 transition between the first 2+ state and the ground state is calculated as B(E2: 
2f of) = 11.8 e^fm^, which agrees with the experimental value of 10.04±1.2 e^fm^. 
Furthermore, the interband transition is calculated: B(E2: 2i — ^ 2^) = 3.99 e^fm"^. If the 
system is axially symmetric, this transition between different i^-values is more suppressed. 
Therefore, the value indicates the 2~^ states have component of a triaxial deformation. 



Using Davydov-Filippov model ||T2[, we can estimate the degree of the triaxiality as a 
Lction of the 7 angle. 
Filippov model as follows: 



function of the 7 angle. The ratios ^(£2- 2+^0+) b(E2' 2+^0+) given in Davydov- 



3-2sin^(37) 



B(E2 : 2^ ^ Ot) _ v/9-8sin^(37) 

B(E2 : 2| ^ 0^) ~ 1 + 3-2.111^(87) 

-y9-8sin2(37) 



+ +\ 20 sin^(37) 

B(E2 : 22 ^ 2i ) _ 7 9-8sm2(37) 



(20) 



B(E2 : 2^ ^ 0^) 1 + 3-2sm^(37) 

-^9-8sin2(37) 

These values are compared with our calculation in Fig. 3. 



(21) 



Fig. 3 

B('E2' 2,'^ ^■2~'~'l 

The solid line in Fig. 3 shows the ratio ) ' + — calculated with the Davidov-Fihppov 
model. The ratio becomes 0.34 in our calculation, and it has a crossing point with Davydov- 
Filippov model around 7 = 19°. The dotted line in Fig. 3 shows the ratio ^(£2- 2+^0+) 
calculated with the Davidov-Filippov model, and it crosses with our result of 0.059 around 
7 = 17° and 22°. These results strongly suggest that ^°Be has a triaxial deformation of 
7 = 15° ~ 20°. Although the a-a core is of axial symmetry and electric charge are only in 
the a's, the recoil effect gives rise to a change from the axial symmetry to the triaxial shape. 
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The 7- value can be deduced also in terms of Q- invariant model. In Ref . , the 7- value 
is related to the reduced matrix elements of Q operator: 



92 = E(otllQII2+)(2+||g||o+), 



(22) 



?3 = 




5:(o+||g||2+)(2+||Q||2;)(2;iiQ||o+), 



(23) 



— 3- = COS 37. 



(24) 



Using these relations and taking the sum over indices i and j up to 2, our calculated results 
correspond to 7 = 20.1°. 

We discuss the electro-magnetic transition between the two 2+ states (i^-mixing effect 
between the two 2"*" states) by artificially changing the strength parameter V^^* in Eq. (15) 
for the spin-orbit term. The Majorana parameter M for the central term is simultaneously 
changed to keep the calculated binding energy constant. In Table III, the B(E2: 2^ 2^) 
values are listed together with the calculated energies and magnetic dipole moments. Using 
the original interaction, the B(E2: 2i 23") value is predicted as 3.99 e^fm'^. 



When we increase the V^'* value, the excitation energy of the second state becomes higher, 
since one of the valence neutrons repulsively feels the spin-orbit interaction. The original 
interaction (V^J" = 2000 MeV and M = 0.6) gives = 5.70 MeV for the 2^ state, but 
the interaction with = 2500 MeV and M = 0.61 gives = 6.82 MeV. The B(E2: 2^ 
value then decreases from 3.99 e^fm^ to 2.20 e^fm'^. Therefore, it is considered that 
with increasing LS strength, the K quantum number of each 2~^ state approaches to a good 
number, where the E2 transition between the two 2"*" states is suppressed. On the other 
hand, when the spin-orbit interaction becomes weaker, the transition rapidly increases. The 
Vq^ value of 1500 MeV gives the B(E2: 2f 2^) values of 9.53 e^fm"^, and when we adopt 
V^'* = 1000 MeV, the value becomes 17.53 e^fm^. Here, the orbits of the valence neutrons 
deviate from ones with good K quantum numbers, and a triaxial a-fa-f-di-neutron clustering 
configuration where the two valence neutrons are strongly correlated becomes important. 

We can intuitively interpret this behavior as follows: when the spin-orbit is weak enough, 
the two valence neutrons form di-neutron, in which the attractive interaction between them 
strongly contributes. However, when the spin-orbit interaction significantly acts and be- 
comes more important, the di-neutron is broken and each valence neutron rotates around 
the core in opposite direction with definite i^T-values. Note that the spin-orbit interaction 
does not act to di-neutron with S = 0. This is close to jj-coupling picture and axial 
symmetry of the system is restored. 

This situation can be interpreted also from nuclear SU3 model. At the SU^ limit (pa- 
rameters a, b, and d 0), $((3/2~)^) (dominant configuration of 2^) and $(3/2^ 1/2") 



Table HI 
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(dominant configuration of 2^) correspond to linear combinations of two SU^ configurations: 
{n^,ny,n,) = (2,0,4), (A, /i) = (2,2) and {n,,,ny,n,) = (1,1,4), (A, ^) = (3,0) (n, = 4 is due 
to the a-clusters along the 2;-axis). This is because, the valence neutrons in 3/2~ and l/2~ 
are expressed as linear combinations of p-orbitals along the x- and y-direction. Here, the 
spin-orbit interaction determines the mixing ratio of the two SU^ representations. When 
it is strong enough, since the jj-coupling picture works well, these two configurations are 
mixed to the almost same amount. However, when the spin-orbit interaction is weakened, 
the binding-energy gain due to the mixing amplitudes changes. Here, the (A, n) = (2, 2) 
configuration, which is a di-neutron configuration, becomes much more important than the 
(A, n) — (3, 1) configuration. This means that a triaxial deformation is induced when the 
spin-orbit interaction becomes weaker. 

Next, we examine the magnetic dipole moment as a probe to determine the strength of 
the spin-orbit interaction. In Table 111, the magnetic dipole moments of these two 2+ states 
are listed, and we predict jj, = 0.72/xjv and OAS/in for the first and the second 2+ states, 
respectively. Unfortunately, the dependence of the total magnetic moment of these states 
with respect to the spin-orbit strength is rather monotonic. However, each component has 
a significant dependence. In Table 111, proton-orbital part and neutron-spin part are also 
listed (neutron-orbital part is of course zero, and proton-spin part also becomes zero due to 
the assumption of the a-clusters) . 

As for the proton-orbital part, using the original interaction (Vq^ — 2000 MeV, M — 0.6), 
2f is calculated to have a larger value than 22 the 2^ state has 1.13 /ij^ and the 22 state has 
0.46 /iAf. The 2f state is mainly of i^' = 0, and therefore the rotation of the a-a core is the 
main source of the angular momentum ( J = 2). However, in the case of 22 , since the valence 
neutrons already have K = 2, the rotation of a-a is not necessary to construct J = 2, and 
hence the proton-orbital part is smaller than 2^. These values depend on the strength of the 
spin-orbit interaction. With decreasing spin-orbit strength, the difference of proton-orbital 
part between these two states becomes smaller due to the increase of the fT-mixing effect 
between the two 2+ configurations, which we discussed. When the strength of the spin-orbit 
interaction is Vq'^ = 1500 MeV, 2i has 1.04 and 22 has 0.65 iin- As for the neutron-spin 
part, in Vq^ — 2000 MeV case (original interaction), 2'^ has —0.41 //jv and 22 has —0.04 //jv. 
This result suggests that the 2i state has both spin-singlet and spin-triplet components of 
the valence neutrons. The 22 state is dominated by the spin-singlet component. In 22 , one 
of the valence neutrons repulsively feels the spin-orbit interaction, and when the two valence 
neutrons construct spin- singlet, it helps to reduce the contribution of the repulsive spin-orbit 
interaction. This is the main reason for the smaller neutron-spin part in 22 . On the other 
hand, if we use weaker spin-orbit strength, the difference of neutron-spin part between the 
two 2"*" states becomes much smaller. When the strength of the spin-orbit interaction is 
V^^' = 1000 MeV, 2t has -0.13 and 2^ has -0.12 /i^v. 

Finally, we comment that the K-mixing effect is also important for the 4"'" states. The 
4+ state with X = at -47.6 MeV in Fig. 2 becomes -50.0 MeV after the X-mixing, and 
the 4+ state with X = 2 at -46.1 MeV is pushed up to -43.5 MeV. The first 4+ state has 
the squared overlap of 0.78 with the K — while 0.63 with the K = 2 state (here, 7^ = 
and K = 2 are not orthogonal). Therefore, the first 4+ state has almost equal contribution 
of the K = Q and the K = 2 basis states. This is because, in the case of 4+, the angular 
momentum vector with K — Q and K — 2 are not necessary to be spatially orthogonal. 
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different from the case of 2"*". Therefore, the i^-mixing for the 4"*" state of ^°Be is very 
strong, and the electro-magnetic transition probabihty among the yrast band deviates from 
a simple rigid-body picture. The B(E2: if 2i) value is calculated to be 11.1 e^fm"^, even 

smaller than the B(E2: 2^ 0^) value of 11.8 e^fm^ and the l^^l t\^S\ I'atio of 0.94 
much deviates from the rigid-body limit of 1.43. If we restrict the model space to = 0, 
then the B(E2: Af — > 2^^) becomes 14.7 e^fm^, and the ratio ^(£2- 2+^0+) '^^ ^'^^ becomes 
closer to 1.43. In nuclear SU3 model, the 4+ states with K = ((A,/i) = (2,2)) and K = 2 
((A,/i) = (3,1)) become the same representation. This character partially remains as a 
strong i^'-mixing effect in the present MO model. 



IV. SUMMARY AND CONCLUSION 

We have applied the a+a+n+n model to ^°Be and discussed the triaxial deformation 
of this nucleus. The orbits for the valence neutrons have been introduced based on the 
molecular orbit (MO) model. In the present model, the spatially extended motion of the 
valence neutrons around the a-clusters are described by linear combinations of Gaussians, 
and the centers of the Gaussians are variationally determined. 

The calculated energy levels show the appearance of two band structures, which have 
dominantly K = and K = 2 components, when the two valence neutrons occupy K'" = 
3/2~ or if = 1/2^ orbits. The first 2+ state has the squared overlap with the K = state 
by 0.96, and the second 2+ state has the squared overlap with the K = 2 state by 0.92. Since 
the second 2"*" state has the component of ii" = by 8%, the electro-magnetic transition 
from the 2^ state (mainly = 0) to the 22 state (mainly = 2) is allowed (3.99 e^fm^). 
Using Davydov-Filippov model |T^, we estimated the triaxiality for the 2+ state. The ratio 

calculated with the Davidov-Filippov model and our model coincide around 

7 = 19°, and B(E2) ratio ^(£2- 2+^0+) indicates 7 = 17° ~ 22°. Originally, the a-a core 
is introduced to be of axial symmetry. However, because of the recoil effect of the valence 
neutrons, the charge distribution deviates from the axial symmetry, and the system becomes 
triaxial. We also discussed the triaxial deformation by artificially weakening the spin-orbit 
interaction. With decreasing spin-orbit interaction, the orbits of the valence neutrons deviate 
from the jj-coupling limit, and the di- neutron configuration becomes important. Here, the 
system becomes 3 body-like and the B(E2 : 22 —>■ 2f) value drastically increases. 

Similar di-neutron component is discussed in weakly bound systems with the so-called 
halo structure. For example, in ^He, in addition to the shell model-like space, the model 
space of di-neutron+^He has been shown to be important [^. This means that a locally- 
correlated di-neutron wave function is important for the description of the valence neutrons 
with small binding energy and spatially extended distribution. It is consistent with our 
discussion on the effect of varying the strength of the spin-orbit interaction. Namely, when 
the valence neutrons with low-binding energies have halo structure, the contribution of the 
spin-orbit interaction between the core and the valence neutrons becomes weak. Here, the 
valence neutrons construct a di-neutron pair with spin singlet , by which they can increase the 
spatial overlap between them and the contribution of the attractive interaction. Therefore, it 
is very challenging to explore the triaxial deformations in ^^Be and ^^Be, which have weakly 
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bound neutrons and also deformed cores. 
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FIGURES 



FIG. 1. The intrinsic density of $((3/2 )^) ((a): xz-plane, (b): xy-plane), where the a-a 
distance is chosen to be optimal 3 fm. 

FIG. 2. The calculated energy levels of ^"^Be. The left and the middle column show the levels 
calculated with the constraints of X = and K = 2. The right column shows the calculated energy 
levels after i^T-mixing. The calculated a+a+n+n threshold energy (—55.0 MeV) is shown as the 
dashed line. Experimentally, the 2^ and the 23" states are observed at = 3.3 MeV and = 5.9 
MeV, respectively. 

FIG. 3. B(E2) ratios p^fo o+^n+! (solid line) and "t+^^+l (dotted line), as a function of 

7 (degree). Our results of 0.34 and 0.059 cross with Davydov-Filippov model around 15" ~ 20°. 
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TABLES 



TABLE I. The values of a, b and the intrinsic energy as a function of the a-a distance. Param- 
eter a for the spin-up and the spin-down valence neutrons ai and 02, respectively, and parameter b 
for them is bi and 62- The wave function of ^°Be is $((3/2-)2) (upper panel), ^{{1/2- f) (middle 
panel), and $(3/2~, l/2~) (lower panel). 



a-a (fm) 


ai = 02 (fm) 


b (fm) 


intrinsic (MeV) 


2 


1.09 


0.99 


-52.09 


3 


1.39 


1.08 


-50.91 


4 


1.71 


1.35 


-44.64 


5 


2.12 


1.53 


-38.04 



Q-n (fm) 


(li = 0/2 (fm) 


b (fm) 


inlrinsic (MoV) 


2 


1.33 


2.05 


-38.06 


3 


1.49 


2.11 


-39.30 


4 


1.63 


2.21 


-35.42 


5 


1.89 


2.35 


-30.34 



a-a (fm) 


ai (fm) 


a2 (fm) 


61 (fm) 


62 (fm) 


intrinsic (MeV) 


2 


1.12 


1.18 


0.83 


1.87 


-50.20 


3 


1.33 


1.37 


1.12 


1.94 


-49.81 


4 


1.66 


1.59 


1.40 


2.09 


-44.16 


5 


2.05 


2.36 


1.59 


2.23 


-37.79 



TABLE II. The electro-magnetic transition probability (B(E2)) in ^''Be after the X-mixing. 
All units are e^fm^. 



B(E2: 




-0+) 


11.8 (Exp. 10.04±1.2) 


B(E2: 


2t- 




0.70 


B(E2: 




-2+) 


3.99 


B(E2: 


4- 


-2+) 


11.1 



14 



TABLE III. The electro-magnetic transition probability (B(E2)) from first 2"*" state to second 
2"'" calculated by changing the strength of the spin-orbit interaction (Vg*). The Majorana parameter 
(M) for the central interaction is also changed to keep the calculated binding energy constant. The 
original interaction is Vq* = 2000 MeV and M = 0.6. Magnetic dipole moment is also predicted. 
nipl) and fi{ns) represent proton-orbital part and neutron-spin part, respectively. 



Fj" (McV) 


1000 


1500 


2000 


2500 


M 


0.58 


0.59 


0.60 


0.61 


0^ (MeV) 


-60.48 


-60.38 


-60.51 


-60.75 


2+ (MeV) 


-57.69 


-57.30 


-57.26 


-57.33 


2^ (McV) 


-56.71 


-55.78 


-54.88 


-53.93 


B(E2: 2^ ^ 2^) (e^fm^) 


17.53 


9.53 


3.99 


2.20 


(2|) M 


0.57 


0.73 


0.72 


0.69 


inipl), ^(ns)) 


(0.70, -0.13) 


(1.04, -0.30) 


(1.13, -0.41) 


(1.16, -0.48) 


// (2+) M 


0.91 


0.59 


0.48 


0.43 


inipl), n{ns)) 


(1.03, -0.12) 


(0.65, -0.06) 


(0.51, -0.04) 


(0.46, -0.03) 
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(a) 




X-ax i s (f m) 



(b) 

3 




X-ax i s (f m) 




K = 



K = 2 



K-m i X i ng 




5 10 15 20 25 30 

Y value (degree) 



